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Abstract 

We study the linearized stability of n- vortex (n G Z) solutions of the magnetic Ginzburg-Landau 
(or Abelian Higgs) equations. We prove that the fundamental vortices (n = ±1) are stable for all 
values of the coupling constant, A, and we prove that the higher-degree vortices (|n| > 2) are stable 
for A < 1, and unstable for A > 1. This resolves a long-standing conjecture (see, eg, [|JT[ ). 
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1 Introduction 



In this paper, we determine the stabihty of magnetic (or Abehan Higgs) vortices. These are certain 
critical points of the energy functional 

E{^,A) = \ j^^ jlVAV'P + (V X Af + ^(l^p - (1) 

for the fields 

A .B? ^B? and ^ : ^ C. 

Here = V — is the covariant gradient, and A > is a coupling constant. For a vector. A, 
V X A is the scalar 9i^2 ~ ^2^1^ and for a scalar V x ^ is the vector {— 82^,81^). Critical points 
of E{ip, A) satisfy the Ginzhurg- Landau (GL) equations 

-AaV + ^(|^|'-i)V' = o (2) 

V X V X A - 9(V^VA'i/') = (3) 

where = Va • Va- 

Physically, the functional A) gives the difference in free energy between the superconducting 
and normal states near the transition temperature in the Ginzburg-Landau theory. A is the vector 
potential (V x ^ is the induced magnetic field), and tp is an order parameter. The modulus of "0 is 
interpreted as describing the local density of superconducting Cooper pairs of electrons. 

The functional A) also gives the energy of a static configuration in the Yang-Mills-Higgs 

classical gauge theory on R^, with abelian gauge group U{1). In this case ^ is a connection on the 
principal U{1)- bundle R^ x C/(l), and il) is the Higgs field (see [lIT|| for details). 

A central feature of the functional E{ip, A) (and the GL equations) is its infinite-dimensional 
symmetry group. Specifically, -E(i/;, A) is invariant under U{1) gauge transformations, 

^ e^^V (4) 

A^ A + V-f (5) 
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for any smooth 7 : R. In addition, E{tp, A) is invariant under coordinate translations, and 

under the coordinate rotation transformation 

V'(x) I— > 'il)[g~^x) ■^{^) ^ 9^{9~^^) (6) 

for g G SO {2). 

Finite energy field configurations satisfy 

IV'I — > 1 as l^l — > 00 (7) 

which leads to the definition of the topological degree, deg(^), of such a configuration: 

: Si ^ S^^ 



deg(V') = deg — 



\x\=R / 

{R sufficiently large). The degree is related to the phenomenon of flux quantization. Indeed, an 
application of Stokes' theorem shows that a finite-energy configuration satisfies 



deg(V) = ^ I (V X 
27r 7r2 



We study, in particular, "radially-symmetric" or "equivariant" fields of the form 

^^^){x) = fn{r)e'''^ = n^^x^ (8) 

r 

where {r,9) are polar coordinates on R^, x-*- = ^(— X2,a;i)*, n is an integer, and 

fn,an : [0,00) R. 

It is easily checked that such configurations (if they satisfy (|^)) have degree n. The existence of critical 
points of this form is well-known (see section |2.1| ). They are called n-vortices. 
Our main results concern the stability of these n-vortex solutions. Let 

= Hess S(^("),A(")) 
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be the linearized operator for GL around the n- vortex, acting on the space 

X = L^{B?, C) © l2(R2, r2). 
The symmetry group of E{iJj^A) gives rise to an infinite-dimensional subspace of ker(L^")) C X (see 



section 3.2), which we denote here by Zsym- We say the n- vortex is (linearly) stable if for some c > 0, 

sym 

and unstable if L*^"^ has a negative eigenvalue. The basic result of this paper is the following linearized 
stability statement: 

Theorem 1 1. (Stability of fundamental vortices) 
For all A > 0, the zizl-vortex is stable. 

2. (Stability /instability of higher- degree vortices) 
For \n\ > 2, the n-vortex is 

stable for A < 1 
unstable for A > 1. 

Theorem |l] is the basic ingredient in a proof of the nonlinear dynamical stability /instability of 
the n-vortex for certain dynamical versions of the GL equations. These include the GL gradient flow 
equations, the Abelian Higgs (Lorentz-invariant) equations, and the Maxwell equations coupled to 
a nonlinear Schrodinger equation. These dynamical stability results are established in a companion 
paper (p|). 



The statement of theorem || was conjectured in |JT| on the basis of numerical observations 
(see |J^). Bogomolnyi (g]) gave an argument for instability of vortices for A > 1, |n| > 2. Our result 
rigorously establishes this property. 

The solutions of (§-§) are well-understood in the case of critical coupling, A = 1. In this case, the 
Bogomolnyi method ( Q ) gives a pair of first-order equations whose solutions are global minimizers of 
E{il), A) among fields of fixed degree (and hence solutions of the the GL equations). Taubes (p?l|, |T2[ ) 



has shown that all solutions of GL with A = 1 are solutions of these first-order equations, and that for 
a given degree n, the gauge-inequivalent solutions form a 2|n|-parameter family. The 2\n\ parameters 
describe the locations of the zeros of the scalar field. This is discussed in more detail in |JT| (see also 
[ BGP[ ) and section ^. We remark that for A = 1, an n- vortex solution (^) corresponds to the case 
when all \n\ zeros of the scalar field lie at the origin. 

The remainder of this paper is organized as follows. In section § we describe in detail various 
properties of the n-vortex. In particular, we establish an important estimate on the n-vortex profiles 
which differentiates between the cases A < 1 and A > 1. In section |^, we introduce the linearized 
operator, fix the gauge on the space of perturbations, and identify the zero-modes due to symmetry- 
breaking. Sections ^ through ^ comprise a proof of theorem |l]. A block-decomposition for the linearized 
operator is described in section ^. This approach is similar to that used to study the stability of non- 
magnetic vortices in [ 0S1[| and |G1]. In section |5|, we establish the positivity of certain blocks (those 
corresponding to the radially-symmetric variational problem, and those containing the translational 
zero-modes) for all A, which completes the stability proof for the itl-vortices. The basic techniques 
are the characterization of symmetry-breaking in terms of zero-modes of the Hessian (or linearized 
operator), and a Perron- Frobenius type argument, based on a version of the maximum principle for 
systems (proposition |6|), which shows that the translational zero- modes correspond to the bottom of 
the spectrum of the linearized operator. A more careful analysis is needed for \n\ > 2. This requires 
us to review some aspects of the critical case (A = 1) in section |6|. The stability/instability proof for 
|n| > 2 is completed in section |^. We use an extension of Bogomolnyi's instability argument, and 
another application of the Perron-Frobenius theory. 

Acknowledgment: the first author would like to thank the Courant institute for its hospitality during 
part of the preparation of this paper, and especially J. Shatah for some helpful discussions. Part of this 
work is toward fulfillment of the requirements of the first author's PhD at the University of Toronto. 
The second author thanks Yu. N. Ovchinnikov for many fruitful discussions. 
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2 The n-vortex 

In this section we discuss the existence, and properties, of n-vortex solutions. 

2.1 Vortex solutions 

The existence of solutions of (GL) of the form (P) is well-known: 

Theorem 2 (Vortex Existence; ||P|, BC| ) For every integer n, there is a solution 



V;W(rc) = /„(r)e™^ ^("Vx) = n^^x^ (9) 

r 

of the variational equations ^)-^- In particular, the radial functions (fn, cin) minimize the radial 
energy functional 

E^-\f,a) = \ 1^°° |(/0^+n^ ^^";/^' +n^^ + \{f^ - lf\^ rdr (10) 

(which is the full energy functional ^ restricted to fields of the form ^) in the class 

{/, a : [0, oo) ^ R | 1 - / G Hi{rdr), - G Ll^{rdr), - G L^rdr)}. 

r r 

The functions fn, o-n are smooth, and have the following properties (for n ^ 0): 

1. < fn < 1, < On < 1 on (0, oo) 

2. f'n,a'^>Q 

3- fn ~ cr", On ~ dr'^ , as r —>■ (c> and d > are constants) 

4- 1 — /n, 1 — a„— >Oasr^oo, with an exponential rate of decay. 

We call (V'^"), ^(")) an n-vortex (centred at the origin). 

It follows immediately that the functions /„ and a„ satisfy the ODEs 

- A./„ + ^'^^"""^' a + ^{f'n - l)fn = (11) 



and 



-/2(l-aO = 0. 



(12) 



Remark 1 To our knowledge, it is not known if solutions of the form ^ are unique. In the appendix, 
we show that for X > Ir? , any such solution minimizes e'^\ 

Remark 2 The functions fn and an also depend on X, but we suppress this dependence for ease of 
notation. When it will cause no confusion, we will also drop the subscript n. 

Remark 3 The discrete symmetry ^ i-^ •ip, A ^ —A of ( GL) interchanges (V'^"\ A^^^) and {^^~'^\ A^~'^^) 
Thus, we can assume n > 0. 

2.2 An estimate on the vortex profiles 

The following inequality, relating the exponentially decaying quantities /' and I — a, plays a crucial 
role in the stability/instability proof. 

Proposition 1 We have 

' fir) > ^^^^f{r) for A < 1 ^^^^ 
/'(r) < ^^^i^/(r) for A>1 

Proof: Define e(r) = f'{r) — "(^~"(^)) J(r). The properties listed in theorem |2| imply that e(r) ^ as 
r ^ and as r — > cxd. Using the ODEs ((|ll])-(|l^)) we can derive the equation 



(_A, + a)e + ^e' = (l-A)/Y 



where 

, , l+n(l — a), rf, na' 
a{r) = 1 ^ 1 + ^ +/' + —> 

and the result follows from the maximum principle. □ 



3 The linearized operator 

In this section, we introduce the hnearized operator (or Hessian) around the n-vortex, and identify 
its symmetry zero-modes. 

3.1 Definition of the linearized operator 

We work on the real Hilbert space 



X = l2(R^ C) © l2(r2; 



with inner-product 




We define the hnearized operator, L^^^ (= the Hessian of E{^lJ,A)) at a solution {ip,A) of 
through the quadratic form 



dedS 



E{^j; + e^ + dv,A + eB + 6C)U=s=o = ((??, C)L^,a{^, B))x 



for all {i,B), {r],C), G X. The result is 




3.2 Symmetry zero-modes 



We identify the part of the kernel of the operator 



which is due to the symmetry group. 



Proposition 2 We have 
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1. 

for any 7 : R 

2. 




(14) 



\ ' =0 (15) 



for 3 = 1,2. 

Proof: We use the basic result that the generator of a one-parameter group of symmetries of E{ij), A), 
apphed to the n- vortex, hes in the kernel of L^"^ . The vector in (|T^ ) is easily seen to be the generator 
of a one-parameter family of gauge transformations (§-||) applied to the n- vortex. Similarly, the vector 
in ( |T5|) is the generator of coordinate translations applied to the n-vortex. □ 

Remark 4 Applying the generator of the coordinate rotational symmetry ^ to the n-vortex gives us 
nothing new, it is contained in the gauge- symmetry case. 

We define Zsym to be the subspace of X spanned by the zero-modes described in proposition |2|. 
We recall that the n-vortex is called stable if there is a constant c > such that 

L^l^x >c, (16) 

and unstable if L^""^ has a negative eigenvalue. 



sym 



3.3 Gauge fixing 

In order to remove the infinite dimensional kernel of L*^"^ arising from gauge symmetry, we restrict 
the class of perturbations. Specifically, we restrict L(") to the space of those perturbations (^, B) & X 
which are orthogonal to the gauge zero- modes (p^). That is. 




for all 7. Integration by parts gives the gauge condition 



Q{'ljjin)^) =\/ -B. (17) 

As is done in we consider a modified quadratic form L^"), defined by 

< a, U^'^a >=< a, L^'^^a > + j - V • 5)^ 

for a = B) G X. Clearly, L*^") agrees with L^") on the subspace of X specified by the gauge 
condition (p!7[). This modification has the important effect of shifting the essential spectrum away 
from zero (see (26)). A straightforward computation gives the following expression for 



l(n) 



^ i\ ( [-A^ + |(2|V|2-l) + i|Vp]e + i(A-l)V2C"+2iVA^.i? 



\B J y 29[V^^C] + [-A + |V^|2]S 

To establish theorem |l], it suffices to prove that L^") > c > on the subspace of X orthogonal to the 
translational zero-modes (|l5|). 

is a real-linear operator on X. It is convenient to identify L^(R^; R-^) with L^(R^; C) through 
the correspondence 

^ Bi 
B2 ^ 

and then to complexify the space X X = [L^(R^; C)]^ via 



B 



B'' = Bi- iB2 



(18) 



iC,B)^iC,^,B'^,B'^). 
As a result, Z^") is replaced by the complex-linear operator 



(19) 



-in) 



diag {-Aa, -Aa, -a, -A} + 



where 



2(2|V^|2-l) + i|V;|2 U^-lW -i{d*A^) i{dA^) 



i{d\i,) 

-i{dAi^) 



i{dA^p) 
-iid*A^) 







'A' 
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Here we have used the notation 

dA = dz- iA 

where dz = di — 182 (and the superscript c has been dropped from the complex function A obtained 
from the vector- field A via (18))- 

The components of F*^"^ are bounded, and it follows from standard results (| RSII[ ) that L is a 
self- adjoint operator on X, with domain 



[^2(R=^;C)]^ 



4 Block decomposition 

We write functions on in polar coordinates. Precisely, 

X = [L2(R2; C)]^ = [Ll, L\S'- C)]' (20) 

where L^^^ = L^(R+,rdr). 

Let pn '■ U{1) Aut{[L'^ {S^ ; C)]^) be the representation whose action is given by 

Pn{e''){tv,B,C){x) = (e^'^^^,e-*"^,e-*^S,e*^C)(i?_ex) 

where Ra is a counter-clockwise rotation in R^ through the angle a. It is easily checked that the 

-in) ~(n) 

linearized operator L commutes with pn{g) for any g G U{1). It follows that L leaves invariant the 

~(n) 

eigenspaces of dpn{s) for any s E iR = Lie{U{l)). The resulting block decomposition of L , which 
is described in this section, is essential to our analysis. In particular, the translational zero-modes 
each lie within a single subspace of this decomposition. 

4.1 The decomposition of L^^'' 

In what follows, we define, for convenience, 6(r) = IiiL_JJ^Il)) _ 
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Proposition 3 There is an orthogonal decomposition 



meZ 



z{n) 

under which the linearized operator around the vortex, L , decomposes as 



where 



with 



and 



LW = -A,(/d) + yir) 

Vj^) = i- diag {[m + n(l - a)f , [m - n(l - a)]^ [m - 1]^, [m + 1]^} + V' 



v 



I 1(2/2 _ 1) + 1/2 1 _ i)/2 // _ ^/ + ^^ 

/'-f*/ -[/' + &/] f 



(21) 



(22) 



Proof: The decomposition (|2l| ) of X follows from the usual Fourier decomposition of L (S ; C), and 

r(n) 

the relation (pO|). An easy computation shows that L preserves the space of vectors of the form 



i{m+n)9 ^Mm—n 
1 '/^ 



-lae 



iim-l)e ^^p^i{m+i)e 



(23) 



and that it acts on such vectors via (|22D. □ 

It follows that Lm^ is self-adjoint on [L'^adl^- ^^^^ ^^^^ convenient to work with a rotated 



version of the operator Lm\ 



Li-) 
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where 



V2 



110 
-110 
11 
1-1 



V2 



110 
1-10 
11 
1-1 



We have 



where 



/ ^ + 6^ + 1(3/2-1) 



-2\m\^ 
-2bf 




-2|m|^ 

^ + b^ + Uf-'^)+f 


-2/' 



-26/ 






-2/' 



m^+l _|_ j2 2 

-2^ ^ + /V 



4.2 Properties of L^^^ 

Proposition 4 /laue the following: 
1. 



3. For \n\ = 1 and m> 2, 



with no zero-eigenvalue. 



ae..(LW) = [inin(l,A),oo) 
L(^) - lS") > 



(24) 



(25) 
(26) 
(27) 



Proof: The first statement is obvious. The second statement follows in a standard way from the fact 
that 

lim V^\r) = diag {A, 1, 1, 1} 
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To prove the third statement, we compute 

L(^) - lS"^ = ^^^^^^ diag {m + 2n(l -a), m - 2n(l - a), m - 1, m + 3} 
which is non-negative, with no zero-eigenvalue for m > 2, n = 1. □ 



Remark 5 In light of (pqj, Twe can assume from now on that m > 0. This degeneracy is a result of 
the complexification ^T^) of the space of perturbations. 



4.3 Translational zero-modes 



The gauge fixing (section |3.3| ) has ehminated the zero-modes arising from gauge symmetry. The 
translational zero-modes remain. 

As written in (15), the translational zero-modes fail to satisfy the gauge condition (0). Further, 



they do not lie in L^. A straightforward computation shows that if we adjust the vectors in (|l5|) by 
gauge zero-modes given by ( |T^ ) with j = —Aj, j = 1,2, we obtain 




(VaV)2 

\^ -(V X A)ei 

where ei = (1, 0) and 62 = (0, 1). Ti and T2 satisfy (p^, and are zero-modes of the linearized operator. 
Note also that T±i decay exponentially as |x| 00, and hence lie in L^. 

" (n) 

It is easily checked that Ti ib iT2 lie in the m = ±1 blocks for Lm ■ After rotation by R, we have 

l'-J^It = 



where 



a' a\ 



T = {f',bf,n—,n 
r r 
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5 Stability of the fundamental vortices 

In this section we prove the first part of theorem |l[ Speciffcahy, we show that for some c > 0, L-lj^^^ > c 
for m 1, and L^i^^^\rp± > c. In hght of the discussions in sections 3^, 4^, and 4^, this wiU estabhsh 
the stabihty of the ibl-vortices. 



5.1 Non- negativity of Lq"^ and radial minimization 

Proposition 5 Lq"^ > for all A. 



Proof: 

From the expression (pi) we see that Lq"^ breaks up: 



= iVo e Mo 



(28) 



(abusing notation shghtly) where 



with 



and 



Wo 



Mo = -Ar{Id) + Wo 



-26/ 4 + 



V 



iVn 



-A. + 62 + |(/2_i) + /2 



V 



-2/' 

-2/' -A, + ^ + /2. ^ 

An easy computation shows that Mo is precisely the Hessian of the radial energy, HessE^'^ (see (p!o|)). 

(n) 

Since the n- vortex minimizes Ei- , we have Mn > 0. It remains to show A'^n > 0. We establish the 



stronger result, A^q > 0. Note that 



A^o — GqGq 
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where 

/ 



V 



Or - f'/f f 

f dr + 1/r 

In fact, Go has no zero-eigenvalue. To see this, we first remark that Go is a relatively compact 
perturbation of Go|a=15 due to the exponential decay of the field components. It follows from an 
index-theoretic calculation done in [ |W| , ^, that Go|a=i is Fredholm, with index 0. We conclude that 
the same is true of Gq (for any A). Finally, it is a simple matter to check that Gq has trivial kernel. If 

= 

it follows that 

and hence that /3 = 0, and so ^ = 0. The relation A'^o > follows from this, and the the fact that 

(TessiNo) = [1,00). □ 

5.2 A maximum principle argument 

Removing the equality in proposition ^ requires more work. First, we establish an extension of the 




maximum principle to systems (see, eg, [LM, PA] for related results). We will use this also in the 
proof that the the translational zero- mode is the ground state of (section ^^ ). 

Proposition 6 Let L be a self-adjoint operator on L^(R"';R'^) of the form 

L = -A{Id) + V 

where V is a d x d matrix-multiplication operator with smooth entries. Suppose that L > and that 
for i 7^ j, Vij{x) < for all x. Further, suppose V is irreducible in the sense that for any splitting of 
the set {1, . . . , d} into disjoint sets Si and S2, there is an i ^ Si and a j ^ S2 with Vij{x) < for all 
X. Finally, suppose that = rj £ with r/ > component-wise, and ^ ^ 0. Then either 
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1. C > or 

2. rj = and ^ < 0. 



Proof: We write C = C"*" ~ C with > component- wise, and compute 

Since and have disjoint support, we have 

r.h.s = J2< ^J^^jkCk > - <C,V> < 0. 

Thus we have 

1. = < <e",^r > 

2. 0= <e7,l/,fee^ > for alii/ /c 

Since L > 0, the first of these implies = and hence = r/. So if ^ 0, then ^+ ^ 0. If = 
and ^+ = 0, replace with — in what follows. An application of the strong maximum principle (eg. 



[GT], Thm. 8.19) to each component of the equation 

now allows us to conclude that for each k, either > or = 0. We know that for some k, > 0. 
Looking back at the second listed equation above, and using the irreducibility of V, we then see that 
= for all j. Finally, we can easily rule out the possibility = for some k, by looking back at 
the equation satisfied by ^fc- Thus we have ^ > 0. □ 

5.3 Positivity of 4""^ 

Now we apply proposition |6| to show Mq > 0. The trick here is to find a function ^ which satisfies 
Mq^ > 0. This allows us to rule out the existence of a zero-eigenvector, which would be positive 
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by proposition ^. To obtain such a ^, we differentiate tlie vortex witli respect to tlie parameter A. 
Specifically, differentiation of the Ginzburg-Landau equations with respect to A results in 



MoC = V (29) 



where 




and 







> 0. 



-(n) 



We can now establish 
Proposition 7 For all X, Lq'' > c > 0. 
Proof: We have already shown in the proof of proposition ^, that A^o > and Mq > 0. Hence, due 



to (^81) and (Pq), it suffices to show that NuII{Mq) = {0}. Suppose MqC = 0, C ^ 0. Proposition 



then implies C > (oi' slse take — C)- Now 

0= <MoC,C> = <C,Moe> = <C,ri> > 

gives a contradiction. □ 

Remark 6 Proposition applied to equation (2^ ) also gives ^ > 0. That is, the vortex profit 



ies 



increase monotonically with A. This can be used to show that the rescaled vortex {fn{r/^/\), an{r / ^/X)) 
converges as X ^ oo to (/*,0), where f* is the (profile of) the n-vortex solution of the ordinary GL 



equation: -Arf*+n'^f*/r'^ + {f*'^-l)f* = 0. This result was established by different means in ^ABQ] 



5.4 Positivity of lJ^^^ 



Proposition 8 L^^^ > with non- degenerate zero-eigenvalue given by T. 
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Proof: Let /U = infspecL^^^^ < 0, which is an eigenvalue by (26). Suppose L^^^^S* = /u5. 
Applying proposition ^ to L^^^ — /i (note that Vi satisfies the irreducibility requirement) gives S > 
(or S < 0). Further, fi is non-degenerate, as if ^ were degenerate, we would have two strictly positive 
eigenfunctions which are orthogonal, an impossibility. Now if /i < 0, we have < 5, T >= 0, which is 
also impossible. Thus is a multiple of T, and /x = 0. □ 

5.5 Completion of stability proof for n = ±1 

We are now in a position to complete the proof of the first statement of theorem |l]. By proposition |^, 
^o"^^^ > c > 0. By proposition | and (||), L^^^^It^ > c > 0. Finally, by (H), L^^^ > c' > for 
|m| > 2. It follows from proposition ^ that L^") > c > on the subspace of X orthogonal to the 
translational zero-modes. By the discussion of section 3^, this gives theorem || for n = ±1. □ 



6 The critical case, A = 1 



In order to prove the remainder of theorem ||, we exploit some results from the A = 1 case. 

6.1 The first-order equations 

Following [^], we use an integration by parts to rewrite the energy (|l|) as 

Eii;, A) = \f + [V X A + i(|^|2 - 1)]2 + i(A - l)(|Vp - 1)^} + 7rdeg(V^) (30) 

(recall, since we work in dimension two, V x j4 is a scalar) where deg(V') is the topological degree of 
V', defined in the introduction. We assume, without loss of generality, that deg{;^^ > 0. Clearly, when 
A = 1, a solution of the first-order equations 

= (31) 
Vx^ + i(|^|2-l) = (32) 
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minimizes the energy within a fixed topological sector, deg(V') = n, and hence is stable. Note that we 
have identified the vector- field A with a complex field as in 

The n- vortices @ are solutions of these equations (when A = 1). Specifically, 

n- = kl-f) (33) 
r 2 



and 



/' = nil^. (34) 



In fact, it is shown in |T2| that for A = 1, any solution of the variational equations solves the first- 
order equations (^l|-|32|). 

Beginning from expression (30) for the energy, the variational equations (previously written as (^ 
^)) can be written as 

^*A^A^P] + V[V X A + idV'p - 1)] + i(A - 1)(|^|2 - 1)^ = (35) 

^[aAV'] - id-z[V xA + ^(IV'P - 1)] = (36) 
(here d\ = —dz + lA is the adjoint of Oa)- 



6.2 First-order linearized operator 

We show that the linearized operator at A = 1 is the square of the linearized operator for the first-order 
equations. 

Linearizing the first-order equations ( |3l[]3^ ) about a solution, (^/i, A) (of the first-order equations) 
results in the following equations for the perturbation, q = {(^,B): 

OaC - iBip = 
V X B + 5R(VS0 = 0. 
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Now using —idzB = V x B — i(V ■ B), and adding in the gauge condition (p!7|), we can rewrite this as 



Lict = 



(37) 



where 



' Oa -ill} 
^ -0 -idz 



If we hnearize the full (second order) variational equations (in the form (35-^)) around {iJj^A), 
we obtain 

d\[dAi - iB^] + iB[dA^] + V[V X 5 + 
+C[V xA + i(|V^p - 1)] + i(A - 1)[(|V;|2 - 1)C + 2V'5R(V^O] = 

and 



i^[dAi - iB^j] + i^idAi^] - id^iV xB + 5R('0^)] = 0. 



Proposition 9 When A = 1, these linearized equations can also be written 

L\Lia = 

Proof: This is a simple computation using the fact that the first-order equations (^ll-|32[) hold. □ 
This relation holds also on the level of the blocks. A straightforward computation gives 

]r(n)| p* p 



where 



m 
r 



^ dr-b f f \ 

dr-b -f 

-/ dr + l/r f 

V / ^ dr + l/r J 



22 



6.3 Zero-modes for A = 1 



It was predicted in |^ (and proved rigorously in IQ) that for A = 1, the Unearized operator around 
any degree-n solution of the first-order equations has a 2|n|-dimensional kernel (modulo gauge trans- 
formations). This kernel arises because the Taubes solutions form a 2|n| -parameter family, and all 
have the same energy. The zero-eigenvalues are identified in [Q, and we describe them here. Let Xm 
be the unique solution of 



m 







on (0, oo) with 



and 



for m = 1,2, 



where 



Xm ~ ?^ as r ^ 



Xm ^0 as r — > oo 



, n. Then it is easy to check that 

F±mWm = 



(38) 



fXm 
fXm 
ix'm + mXm/r) 
V -ix'm + ^Xm/r) J 



We remark that 



1 - a 



Xi 



and it is easily verified that for A = 1, Wj-i = T are the translational zero-modes. 



7 The (in) stability proof for |n| > 2 

Here we complete the proof of theorem |^. 
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The idea is to decompose Lm^ into a sum of two terms, each of which has the same (translational) 
zero-mode (for m = 1) as . One term is manifestly positive, and the other satisfies restrictions of 
Perron-Frobenius theory. 

We begin by modifying F^, and defining, for any A, 



/ 



f ' 



fq 







dr-^ -f 

-f dr + 1/r 




/ 





m 
r 



dr + l/r J 



where we have defined 



n(l-a)/ 
rf 



(39) 



W„ 



and dr • q denotes an operator composition. By (34), we have q = 1 for A = 1. We also set, for 
m = 1, . . . , n, 

^ q-'fx^ ^ 

fXm 

Now Wm has the following properties: 

1. W±i is the translational zero- mode T for all A 

2. when A = 1, Wm = Wm, m = ±1, . . . , in, give the 2n zero-modes ( p8[ ) of the linearized operator. 

These Wm were chosen in [Q as candidates for directions of energy decrease (for \m\ > 2) when A > 1. 
Intuitively, we think of Wm as a perturbation that tends to break the n-vortex into separate vortices 
of lower degree. 

Now, Fm was designed to have the following properties: 

1. Fm = Fm when A = 1 (this is clear) 
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2. FmWm = for all m and A (this is easily checked). 
A straightforward computation gives 

l'^^ = F:^Fm + JMrn (40) 

where J = diagjl, 0, 0, 0} and 

Mm = lm- qlmQ + (A - q'^)f^ 

with 



2 



By construction, when m = 1, the second term in the decomposition ( |40| ) must have a zero-mode 
corresponding to the original translational zero-mode. In fact, one can easily check that Mi/' = 0. 

Proposition 10 For \n\ > 2, Mi has a non- degenerate zero-eigenvalue corresponding to f , and 

Ml > A < 1 
Ml < A > 1 

V 

Proof: We recall inequality (|l3|), which implies that for A < 1, g < 1, and for A > 1, g > 1. The 
operator Mi is of the form 

Ml = (1 - g2)(-A,.) + first order + multiphcation. (41) 

One can show that Mi is bounded from below (resp. above) for A < 1 (resp. A > 1). We stick 
with the case A < 1 for concreteness. Suppose Mirj = firj with /x = inf specMi < 0. Applying the 
maximum principle (eg proposition ^ for d = 1) to (41), we conclude that r/ > 0. If /x < 0, we have 
< Vi f 0) ^ contradiction. Thus /x = 0, and is non-degenerate by a similar argument. □ 
We also have 
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Lemma 1 For m > 2, Mm — Mi is non-negative for A < 1, non-positive for A > 1, and has no 
zero- eigenvalue. 

Proof: This follows from the equation 

Completion of the proof of theorem^: Suppose now A < 1. Since Fm^m is manifestly non-negative, 
and Mm > Mi for m > 2, we have Lm^ > for m> 1 (with only the translational 0-mode). Combined 
with ( pBD and propositions ^ and ^, this gives stability of the n- vortex for A < 1. 

Now suppose A > 1. By (^0|), proposition 10 and lemma ||, we have for m = ±2, ...in, 



<Wm,L^m^Wm> < 0. 

We remark that Wm corresponds to an element of the un-complexified space X, and so L^") has 
negative eigenvalues. This establishes the instability of the n-vortex for |n| > 2, A > 1, and completes 
the proof of theorem ||. □ 

8 Appendix: vortex solutions are radial minimizers 

Proposition 11 For A > 2n^, a solution of the equations {T\- li ) minimizes Er ^ . 



Proof: It suffices then to show Mq = HessSl^^ > (see section ^). We write Mq = Lq + Zq where 

I/O = diag{l, — A^} 

with / = -A,. + 6^ + A(/2 _ 

' 2A/2 -26/ 
-26/ 4 + /2 
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We note that If = (one of the GL equations). It follows from the fact that / > and a Perron- 
Frobenius type argument (see 0S1[ ) that / > with no zero-eigenvalue. It suffices to show Zq > 0. 
Clearly tr(Zo) > 0, and 

det(Zo) = 2A/^ + ^[A - 2n'^{l - af] 
is strictly positive for A > 2n^. □ 
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